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We analyze the asymptotic dynamics of quantum systems resulting from large numbers of iter- 
ations of random unitary operations. Although, in general, these quantum operations cannot be 
diagonalized it is shown that their resulting asymptotic dynamics is described by a diagonalizable 
superoperator. We prove that this asymptotic dynamics takes place in a typically low dimensional 
attractor space which is independent of the probability distribution of the unitary operations ap- 
plied. This vector space is spanned by all eigenvectors of the unitary operations involved which 
are associated with eigenvalues of unit modulus. Implications for possible asymptotic dynamics of 
' iterated random unitary operations are presented and exemplified in an example involving random 

controlled-not operations acting on two qubits. 

^ . PACS numbers: 03.65.Ud,03.67.Bg,03.67.-a,03.65.Yz 



O 
O 



0^ 
O 



X 



I. INTRODUCTION 



In recent years the rapid advancement of quantum technology with its capabihties of controlling individual quantum 
systems has given rise to impressive developments in the areas of quantum information science and high-precision 
quantum metrology.^ In particular, current experiments on large ensembles of interacting quantum systems open 
, interesting perspectives to investigate in detail not only the transition from quantum to classical behavior but also to 

■ trace down those quantum phenomena or effects that still are observable on the mesoscopic or macroscopic scale. A 
^ paradigm of such large physical systems are interacting networks whose dynamics is currently investigated intensively 
O^i in the classical domain.^ Such networks are capable of simulating the behavior of real world systems like the internet 

or social dynamicsi^ Typically, in these systems a number of modes representing physical objects are coupled to each 
^-H . other by random interactions. A particularly interesting issue is to determine the dynamics of the system. In view 
^ ' of these current activities the natural question arises which characteristic properties govern the dynamics of such 
networks if each classical node is replaced by a quantum system and, correspondingly, the classical interactions by 

■ quantum operations. 

\ In general, determining the time evolution of large quantum systems is difficult and analytic or closed-form solutions 
. ■ are possible in exceptional cases only. In particular, this applies to the dynamics of open quantum systems in which a 
OO ' large quantum system is in contact with an additional physical system. The influence of such an external system can 
, be taken into account in various ways. In special cases it may be described by randomly applied unitary operations. 
Such a case is realized, for example, if the nodes of a large quantum network represent participants of a quantum 
communication network and if these nodes establish node-to-node communication in a random way by using quantum 
protocols which can be described by unitary transformations. A natural question arising in this context is what is 
the resulting quantum state of the network after a large number of such communication steps. More generally, such a 
• quantum network involving random unitary transformations can characterize the dynamics of any interacting quantum 
' system in which the interactions involved can be described by repeatedly applied random unitary transformations. 
A natural approach to determine the dynamics of a quantum system involves diagonalization of the generator 
of the time evolution. This way the dynamics can be determined in a convenient way even in the asymptotic 
limit of arbitrarily long interaction times. The situation becomes significantly more complicated for open quantum 
systems because the relevant generators are often non-hermitian and not normal*^ so that they cannot be diagonalized. 
Nevertheless, in such cases it is still possible to use the Jordan canonical form (see Appendix |^ of these operators 
for determining the dynamics for arbitrarily long interaction times. This leads to the highly nontrivial problem of 
handling generalized eigenvectors of the relevant generators which are in general not orthogonal. 

Motivated by these aspects in this paper we address the problem of determining general properties of the asymptotic 
dynamics of quantum systems whose dynamics is governed by repeated applications of random unitary transforma- 
tions. A main goal of this paper is to demonstrate that the Jordan canonical form of the generators of random unitary 
transformations have rather unexpected and useful special properties which allow to obtain even closed-form expres- 
sions for the asymptotic quantum state resulting from a large number of iterations of random unitary transformations. 
It will be proved that there is always a vector subspace of so-called attractors on which the resulting superoperator 
governing the iterative time evolution of quantum states can be diagonalized and in which the asymptotic quantum 
dynamics takes place. As a main result a structure theorem is derived for this set of attractors which allows to deter- 
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mine them in a convenient way. Furthermore, it is shown how the asymptotic iterative dynamics of arbitrary quantum 
states can be written in terms of these attractors. Based on these findings we show that in general the asymptotic 
dynamics is non-monotonic. Finally, aspects of these general properties are exemplified by studying in detail the 
dynamics of two qubits which are coupled by randomly applied controlled- not operations. It should be mentioned 
that some of the results characterizing the asymptotic dynamics can also be obtained by a different approach which 
uses special properties of random unitary transformations in order to construct a convenient Ljapunov functioni^ 

This paper is structured as follows. In Section [TT] we summarize basic properties of random unitary transformations 
which are useful for our subsequent discussion. In Section lllll we examine special properties of the Jordan canonical 
form of random unitary maps. The central statement of the paper, namely the structure theorem for attractors of 
random unitary operations, is derived in Section HVl Characteristic properties of attractors are investigated in Section 
fVl Section [VII is devoted to important implications resulting from the structure theorem. Finally, as an example the 
asymptotic dynamics of two qubits which are coupled by random controUcd-not operations is discussed on the basis 
of our general results in Section IVlII 



II. BASIC PROPERTIES OF RANDOM UNITARY OPERATIONS 

A random unitary operation (RUO) $ is a completely positive trace-preserving map admitting a convex decompo- 
sition of the form^ 

m 

Hp) ^Y.P^u,pul (1) 

i=l 

Thereby, Ui denotes a unitary operator acting on a Hilbert space and this unitary operation is applied onto the 
quantum state p with probability > so that ^^^iPi — 1. These latter probabilities take into account classical 
uncertainties in the realizations of the unitary quantum evolution involved. This uncertainty can be the result of 
an unknown error mechanism or of an unknown unitary evolution involving an additional ancillary system. In the 
following we are interested in the asymptotic dynamics resulting from many iterative applications of Starting with 
our quantum system in the initial state p{0), the {n+ l)-st step of this iteration procedure changes the state after the 
n-th iteration p{n) to the state p{n -I- 1) = <f>(p(n)). Our aim is to analyze the asymptotic behaviour of this iteration 
procedure. The random unitary map <& of Eq.([T|) belongs to the class of bistochastic or doubly stochastic maps^i^ 
which leave the maximally mixed state invariant, i.e. 

m 

<f(/)=^KC/,C//=/, (2) 

i=l 

and it acts on the Hilbert space B{Jff) of all linear operators defined on a d-dimensional Hilbert space Jif. The 
dimension of the input and output system is the same. The Hilbert space B{Jif) is equipped with the Hilbert-Schmidt 
inner product {A, B)hs — Tt{A^ B) for all A,B £ I3{H). The adjoint operator of $ is given by 

m 

$t(A)=^K[/Mt/.. (3) 

1=1 

In general, the RUO $ is neither hermitian nor normal and consequently is not diagonalizablc. Therefore, its 
resulting iterated dynamics has to be analyzed with the help of Jordan normal forms'* (see Appendix It is a main 
goal of our subsequent discussion to prove that the Jordan normal forms of RUOs have interesting special properties 
which are particularly useful for the description of their asymptotic iterated dynamics. In particular, there exists a 
Jordan base in the Hilbert space B{J^) in which the matrix of the map ([1]) has a block diagonal form (|Aip . 

Let us formulate first several characteristic properties of RUOs. 

Lemma II. 1 The random unitary map $ defined by the relation (QJj fulfills the following properties: 

1) The norm of the RUO ^ induced by the Hilbert-Schmidt norm of the Hilbert space B{J^) equals unity. 

2) If X is an eigenvalue of the map $, then |A| < 1. 

3) Let X\ e B{J^) be a generalized eigenvector corresponding to the eigenvalue A of the map then X — 1 or 
TiXx = 0. 
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Proofs. 

(1-2) First we prove that the Hilbert-Schmidt norm is unitarily invariant. For this purpose consider an arbitrary 
operator A e B(J^) and two unitary operators U,V ^ B{Jf). As a trace of matrix products is invariant under cycHc 
permutations we get 



\\UAV\\hs = {Tr [{UAVyiUAV)]y - {Tr(AU)}^ = \\A\\hs. 



(4) 



Therefore one can show that ||$|| = sup||^||^g<i ||$(A)||/f5 = 1. Let A e B{Jf ), then the Hilbert-Schmidt norm of 
the operator ^{A) is bounded by 



\mA)\\Hs = 



Y,P,U,AUI 



<Y,pA\u,AU} 



HS 



HS 



\\A\\ 



HS 



(5) 



Moreover, we have ||$(/)||_f/s = H-^llffS- Hence, ||$|| — 1 and consequently |A| < 1. 

(3) If Xx is a generalized eigenvector corresponding to an eigenvalue A of the map <&, then there is a rt g N such that 
((f> — \I)^{X\) = because a simple calculation yields 



Tr{($- A/)"(Xa)} = (1 - A)"TrXA = 0. 
This equation can be fulfilled only if A = 1 or Tr^A = 0. 



(6) 



Thus, all Jordan blocks in the Jordan normal decomposition of the map $ correspond to eigenvalues A with |A| < 1. 
For our subsequent discussion let us introduce the following notation. Suppose that A is an eigenvalue of the map $. 
We denote the corresponding eigen-subspace by Ker(<I> — AJ), i.e. 



Ker($ - A/) = {X <E B{M')\^{X) = \X} , 
and the range of the map $ — A/ by Ran(<l> — A/), i.e. 

Ran($ - A/) = {X e B{.ytf)\3Y e e(^), X = - \Y} . 



(7) 



(8) 



Furthermore, let us define dx — dim(Ker(<i> — A/)) and (T|i| as the set of all eigenvalues of the linear map $ satisfying 
|A| = 1. Finally, the vector subspace spanned by all eigenstates corresponding to eigenvalues A with |A| < 1 we call 
the attractor space of the RUO ^ and denote it by Atr(<i>), i.e. 



Atr($) = Ker(<I' - A/). 



(9) 



We call elements of this subspace attractors of the dynamics because, as we will show later, the asymptotic iterated 
dynamics of the RUO is completely determined by these linear operators. 



III. JORDAN CANONICAL FORM OF RANDOM UNITARY OPERATIONS 



In this section we prove that all Jordan blocks corresponding to eigenvalues A with |A| < 1 are one-dimensional. 
In other words, generalized eigenvectors corresponding to eigenvalues |A| = 1 are all eigenvectors. This statement is 
equivalent to the following theorem (for details see Appendix |^ . 

Theorem III.l Let $ : B{J^) — > B{J^) he a random unitary operation defined by (QP and A its eigenvalue satisfying 
|A| — 1, then we have 

Ker($ - XI) n Ran($ - XI) = {0}. (10) 

Proof. We prove this theorem by contradiction. Suppose there is an operator Q ^ A E B{,ytf) and A E Ker($ — 
XI) n Ran($ - XI). This implies $(A) = XA and there is an operator 0^ B E B{J^) such that $(B) = XB + A. By 
induction one can conclude 



$"(S) = A"B + nA"-M 



(11) 
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and consequently 

n||A|| - \\B\\ < ||A"B + nA"-iA| 
Because the resulting inequality 



|$"(B)||<||$in|i3|| = ||i3| 



<^l|i^ll 

has to be fulfilled for arbitrary n eN the only alternative left is that ^ = 0. 



(12) 



(13) 



Let Yj.fe (j E p, k E {1, 2, . . . , dim(Jj)}) (compare with Appendix E)) be the Jordan basis of the RUO <i>. Jj 
is a Jordan block corresponding to an eigenvalue Xj with a basis formed by the generalized eigenvectors Yj^k {k E 



{1,2,..., dim( Jj)}). Let p(0) E B{Jif) be an input density operator. We denote by Pj I the parameters of the unique 
decomposition of the density operator p(0) E B{Jf) into this basis, i.e. 



P dim(Jj) 

j=i k=i 



(14) 



Consider now the density operator p{n) = $"(^(0)) describing the physical system after n iterations and denote its 

decomposition coefficients p4)) into the same basis by Pj"'/}- It is clear that the coefficients l3j"j} corresponding to 
eigenvectors of eigenvalues A^- S cr|i| evolve simply as 

(15) 



/3(")=A,"/3f . 

(We omit the second index k intentionally because in this case all the Jordan blocks are one dimensional.) 

Now we have to analyze the behavior of the remaining coefficients. It is governed by the following theorem which 
quantifies how the remaining coefficients /^j'^ , corresponding to Jordan vectors Yj,k with |Aj| < 1, evolve. 

Theorem III. 2 Let $ : B{.y/f) — > B{J^/f) be a quantum random unitary operation defined by (J^l with its Jordan basis 
Yj,k (j E p, k E {1, 2, . . . , dim( Jj)}J and p(0) E B{J^) be an input density operator. Furthermore, let f3j"j^ be the 
decomposition coefficients of p(n) ~ $"(/3(0)) into this Jordan basis, i.e. 



p dim(Jj) 

pw-E E f^S^^^^- 

3 = 1 k=l 



(16) 



For any eigenvalue Xg (\Xs\ < 1) of the map $ with its corresponding Jordan block Jg and its Jordan chain Ys k 

n) 
■,k 



(k E {l,2,...,dim(Js)}^ the coefficients /3^"2 vanish in the limit of large n 



lim /3("2 

n^+oo ' ''''' 



0, 



for Vfce {l,2,...,dim(J,)}. 



(17) 



Proof. This theorem follows directly from the fact that the Jordan block (Js)" of dimension dim(Js) with 



J. = 



fXs 1 
A. 

V 



1 

XsJ 



(18) 



vanishes in the limit of large numbers of iterations n, i.e. 

lim (J,)" = 0. 

n— >oo 

One can check that the entry {J^)ij {i < j < dimJs) of the upper triangular matrix (Js)^ fulfihs the ineqnahty 



iw 



{j - i) 



(19) 



(20) 
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so that we obtain the relation 



hm (JDy -0, Vi,j e {l,...,dimJ,}. 



(21) 



In view of this theorem the asymptotic dynamics of the state p{0) under iterations of the random unitary operation 
$ is given completely in terms of its attractors. The remaining coefficients of the decomposition of the initial state 
p(0) (jl6p become vanishingly small after sufficiently many iterations of the map. An interesting question which will 
be addressed in the following is how to determine the set of attractors. 



IV. STRUCTURE THEOREM FOR ATTRACTORS 



Let us now study the structure of the attractors, i.e. of all eigenspaces Ker("I> — XI), with A G cr|i|. In the case 
of random unitary operations the following powerful theorem can be proved which allows us to specify the space of 
attractors of the RUO $. In this context it should be also mentioned that for the more general case of arbitrary unital 
quantum operations interesting general results have been derived by Kribs^'^*' recently. 

Theorem IV. 1 Let $ : B{J^f) — > B{.y/f) be a random unitary map $ (QJ) and X G (T|i|. Then the eigenspace 
Ker($ — XI) corresponding to this eigenvalue A is equal to the set 



Dx := {X e B{Ji^)\UrX = XXU^ for i = 1, . . . ,m} . 



(22) 



Proof. The map $ is unital, that is <&(/) = /. Therefore, every X e D\ fulfils $(X) = XX and thus D\ C Ker($-A/). 
To prove the converse, let us consider X e Ker($ — XI). li X = 0, then X e D\. So let us assume that X ^ 0. Using 
the unitary invariance of the Hilbert-Schmidt norm we get 



|X|| = ||AX|| = 



J2p,u,xuI 



4=1 



<J2P^\\U^XU| 



= 11X11. 



Therefore, the inequality (|23p is in fact an equality and can be rewritten in the form 



(23) 



^PiVt,^PiVt = ^p,{v,,v,)^^ 



i=l 



\i=l 



with Vi — UiXUj. Hence we get 



^ 2pipj{vi, Vi) 2 (wj, Wj) 2 = ^PiPj [{vi, Vj) + {v.j,Vi)] = ^ 2piPjRe{vi,Vj) < ^ 2piPj \Re{vi,Vj 

i<j i<j 

< 2p,pj \ {v^,Vj)\<Y^ 2p,pj {v^,v^)^Vj,Vj)^. 



i<j 



(24) 



(25) 



Because the left and right side of the relation are the same, all inequalities are actually equalities. In particular, 
we have 



Re{vi,Vj) = \{vi,Vj)\ = (vi,Vi)-^{vj,Vj)'^ ^0 for all ij e {1,...,to} 



(26) 



which can be fulfilled if and only if Vj = PijVi (for all i,j) with /3y > 0. From the unitary invariance of the Hilbert- 
Schmidt norm 



\X\ 



UiXU. 



A: 



u.xul 



\x\ 



we conclude that (3ij = 1 for all z, j S {1, ...,to} and hence 

UiXUl = U2XUI = ... = u„,xul. 

Finally, using the equality $(X) = XX we obtain U^XU} = XX, i.e. X G Dx- 

As a consequence of this structure theorem IIV.II the following corollary can be proved. 



(27) 



(28) 
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Corollary IV. 2 The random unitary operation $ defined by ([7p fulfills the following properties: 

1) If X is an eigenvalue of the operation $ fulfilling \X\ = 1, then 

Ker($ - XI) _L Ran($ - XI). (29) 

2) If Xi,X2 are two different eigenvalues of the operation $ fulfilling |Ai| = IA2I = 1, then 

Ker($- Ai/) ± Ker($- A2/). (30) 

Proof. First, from the theorem pVl|) follows that if $(X) = XX and |A| = 1, then ^\X) = X*X and $t(xt) = AX^. 
In order to show that the set Ker($ — XI) is orthogonal to the set Ran($ — XI) we have to prove that {K,R) = is 
fulfilled for arbitrary elements K G Ker(<i> — XI) and R e Ran(<i> — XI). Therefore, there is an operator Q g B{J^) 
with R ^ Y^tPiUiQU} - XQ. Hence, using theorem IIV. II we have the orthogonality relation 

{K,R) = Tr{/^ti?} = Y,P^^^{K^U^QU}} - XTy{K^Q} = A^p,Tr{/^tg} - XTy{K^Q} = 0. (31) 

i i 

The second property is a consequence of the identity 

{X,,X2) ^ ^ iX,,^X2)) = ^ {^HX,),X2) = ^{X^,X2) (32) 

A2 A2 A2 

which is valid for Xi e Ker($ — Xil) and X2 G Ker((f> — A2/) and for any mutually different non-zero eigenvalues Ai 
and A2. Therefore, the last equality can be satisfied only if (Xi, X2) = 0. 

■ 

This corollary together with theorem IIII. II has the following important consequence. 

Theorem IV. 3 Let $ : B{Jif) B{Jif) be a quantum random unitary operation defined by (Q]) and p{0) G B{J^) 
be an input density operator, then the asymptotic iterative dynamics of the state p(0) under the evolution map $ is 
given by 

dx 

Poo(n)- A"Tr{p(0)4jXA,. (33) 

•''£'^111 ,i=i 

and satisfies the relation 

lim ||p(n)-poo(n)|| =0 (34) 

n — >OQ 

with p{n) = 4>"(/3(0)) and with the complete set of orthonormal basis elements Xx,i (i £ {1,2, . . . ,dx}) of the space 
Ker($-A/). 

Proof. In order to prove this theorem we have to show that the mutually orthogonal subspaces 

3o = Ker($ - XI) and 3^l = Ran($ - XI) (35) 

are invariant under the map <I> and that they satisfy the relation 3o © = B{.J^). The second claim is a direct 
consequence of corollary IIV. 21 The first claim follows from the fact that all subspaces Ker($ — XI) and Ran($ — XI) 
are invariant under the map that is, $(Ker(<I> - XI)) C Ker($ - XT) and <i>(Ran($ - XI)) C Ran(<I> - XI). Now we 
can choose some orthogonal basis vectors X\^i in the subspaces Ker($ — XI) with |A| = 1 and the Jordan basis Yj^k of 
the map $ restricted to the subspace 3i. These vectors form a basis of the Hilbert space B{Jf). Now we consider a 
decomposition of p{0) into these basis vectors. As was shown in theorem lIII.2l the part corresponding to the subspace 
3i vanishes for n +00 and the dynamics of the state p{0) on the subspace Jq is given by (jlSp . 
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V. BASIC PROPERTIES OF ATTRACTORS 



In this section we discuss some basic properties of RUOs which are useful for obtaining the complete set of attractors. 
A basic property arises straightforwardly from the theorem lIV.il 

Proposition V.l 1) Let X\-^ and Xx^ be attractors of the RUO (QP corresponding to eigenvalue Ai and A2, re- 
spectively, then the product of these attractors X\-^^X\^ is either an attractor corresponding to eigenvalue A1A2 
or it is the zero operator. 

2) Let Xx he an attractor of the RUO (QJ) corresponding to the eigenvalue \, then X^-^ is also an attractor of the 
RUO (Qp corresponding to the eigenvalue A* . 

Proof. This proposition follows from the identities 

UiXxiXx2 = MXxiUiXx2 = ^iMXxiXx2Ui (36) 

and 

U,XI - {XxU})^ = {\U}Xx)^ = A*4c/„ (37) 
which are valid for all i (z rfi :— {1, . . . , m}. 

■ 

Based on our preceding analysis a single step of the asymptotic dynamics is described by the superoperator 

$as.(.)= E ^Tr(xl,(.))XA,., (38) 

Aecr|i| ,1=1 

which fulfils the property 

hm ||<i>"-<fLJI=0. (39) 

n — >CfO 

The superoperator psp is a unital quantum operation. In order to prove this statement, let us define the projector 

dx 

V{.)= J2 Tr(xL(.))x,,„ (40) 

Ag(T|i| ,1=1 

which projects all elements of the vector space B{J^) onto the attractor space Atr($). The structure theorem llV. 1 1 and 
the orthogonality of all elements of the attractor space ensure that [<&,7^] = and [<I'ass,'P] — 0. These commutation 
properties imply that for any integer n the action of the superoperator $"3^ on an arbitrary operator A G B{.y/f) is 
given by 

KssiA) = {V{A)) ^ UIV{A)U}: (41) 

for an arbitrary G rh. Thus, for any integer n the action of the map $",,5 on the Hilbert space B{J^) is a sequence 
of a projection onto the attractor space and a unitary operation. As a consequence it is a completely positive map 
and in view of Eq. (j39p it describes the dynamics of the iterated random unitary operation in the asymptotic limit of 
large numbers n of iterations. 

It is instructive to analyze this property of complete positivity also from another perspective by using the concept 
of dynamical matrices (compare with Appendix |B|.^ In order to obtain the dynamical matrix of the asymptotic map 
$ass we first calculate its matrix elements in an orthonormal basis, i.e. 

{^ass)Z = {m^,\<^>,,,\nl.) = ATr{(|n)(H)Xl,,}Tr{(|m)(/i|)^X,^.} 

Aecr|i| ,1=1 

- E M^A,.)L(^A,.)„,. (42) 

ASctm 1 ,1=1 
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The elements of the dynamical matrix are defined by 

dx 



{D^.JZ = i'^assTP = E M^A.);. (^A.)™„ (43) 

SO that one obtains the relation 

= E E M(^A..L„H(ni]0[(x,.,);jM)H] = E ^^A.^^^A,.- (44) 

Using the identity (jBSp one can rewrite the d? x c?^ dynamical matrix as an operator acting on dx d matrices according 
to 

dx 

D^.ss{A)= E ^^A.^^^I,,, (45) 

Ae<T|i| ,i=l 

where A is an arbitrary d x d matrix. Expressions (|44p and (|45p describe the same dynamical matrix. The first 
relation describes it as a map acting on reshaped vectors of length d^ and the second one as a map acting on d x d 
matrices (for details see Appendix IB|) . Both expressions are useful to determine the properties of attractors. 

According to Eq. (|B9|) the dynamical matrix (|44|) is always hermitian. Due to proposition IV. II this property is 
fulfilled. Furthermore, the partial trace of the dynamical matrix (I44p over each subsystem yields the identity operator. 
For a RUO of the form of Eq. ([T]) both properties lead to the condition 

di 

^Tr{4jXi,, = /. (46) 

i=l 

SO that the dynamical matrix is positive. With the help of equations (US)) and (IB2p we find that this positivity is 
equivalent to the relation 

''a f 

^ XTT[A^Xx^^AX^>0 (47) 

AgCT|i| ,i=l 

which has to be fulfilled for an arbitrary d x d matrix A. In view of theorem lIV.il we can thus conclude that the map 

dx 

Kssi-)= E A"Tr{(.)xljX,,„ (48) 

Aeo-|i| ,i=l 

is a trace-preserving and completely positive unital map for an arbitrary n G 1^. 



VI. DISCUSSION AND IMPLICATIONS 



Let us summarize and comment the results obtained so far for the asymptotic behaviour of a quantum system under 
a RUO. 

First of all, the asymptotic iterative dynamics is determined completely by the attractor set Atr(<l>) of a RUO. The 
Hilbert space B{J^) can be decomposed as B{J^) = Atr(<i>) © (Atr(<l>))"'" with _L denoting the orthogonal complement 
with respect to B{J^). Both mutually orthogonal subspaces, i.e. Atr($) and 3i = (Atr(<i>))^ are invariant under the 
RUO (P) and we proved that the component of any initial quantum state in the subspace 3i vanishes after sufficiently 
large numbers of iterations. Furthermore, we proved that the vector space of attractors Atr(<i>) is spanned by all 
elements X of the set B{Jif) which fulfil the generalized commutation relations UiX = XXUi for all unitary operators 
Ui of the decomposition ([1]) and for all eigenvalue A with |A| = 1 . 

The calculation of the asymptotic iterated dynamics of the random unitary map ([1]) can be divided into four steps: 
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• One determines the set a\i\. Usually, this step is highly nontrivial and depends significantly on the particular 
unitary Kraus operators involved. Any additional properties concerning the structure of the unitary operators 
Ui involved, for example, simplify this task considerably. In particular, the exploitation of symmetries may be 
useful in this respect. 

• One identifies the set of attractors of the RUO Atr($). This step involves the calculation of all eigenspaces using 
the generalized commutation relations Ker(<i> — XI) — D\ for all A e cr|i|. 

• One chooses an orthonormal basis X\ .i in each subspace Ker($ — XI) for A G cr|i|. 

• One calculates the asymptotic iterated dynamics according to the relation 

p(n) = a>"(p(0))(n»l)= ^ A"Tr{p(0)4jX,, (49) 

-^e(T|i| ,j=i 

which is valid asymptotically for n S> 1. 

These general features imply some important consequences. Firstly, the set of attractors Atr($) and its correspond- 
ing spectrum is independent of the nonzero probabilities pi defining the convex decomposition of the RUO in equation 
([1]). Thus, two RUOs with the same unitary operators in their convex decompositions H]) have the same attrac- 
tors space Atr(<l>). The nonzero probabilities pi determine only how fast an input state converges to the asymptotic 
attractor space. 

Another simple consequence arises if the ensemble of random unitary operators defining the RUO $ contains the 
identity operator / (apart from a global phase). Theorem IIV. II implies that the only possible eigenvalue of the map 
$ is A = 1 . Hence from the set of attractors only fixed points can be formed and the resulting asymptotic dynamics 
is stationary. Moreover, assume that the unitary operators Ui are generators of a finite multiplicative group. As 
any group contains a unit element all possible eigenvalues of the RUO <i> fulfil the relation A"^ = 1 for some integer 
n\ G N. As a consequence the resulting asymptotic dynamics is periodic. Such a periodic asymptotic dynamics is also 
obtained obtained if the unitary operators Ui form an irreducible set of operators, i.e. they have no common nontrivial 
invariant subspace. This can be proven as follows. Consider an eigenvalue A of the random unitary operation ([1]) with 
|A| = 1 and its corresponding eigenvector Xx ^ 0. Using theorem IIV. II it can be checked that Ui{K&x(X\)) C Ker(X>,) 
and U i{R2iX\(X\)) C Ran(XA) is fulfilled for all i G m. Thus, X\ 7^ is an invertible operator. Let a 7^ be an 
eigenvalue of the operator X\ and Ua its corresponding eigenvector. From the equation 

X\UiUa = XUiXxUa ~ aXUiUa (50) 

follows that also aA is the eigenvalue of X\ and UiUa are its corresponding eigenvectors. Therefore also aX^, aX^, ... 
are eigenvalues of X\. Eigenvectors corresponding to different eigenvalues are linearly independent. Therefore there 
is n g N such that A" — 1. Moreover, the direct sum of all eigensubspaces corresponding to eigenvalues X^a 
(j € {0, 1, ...,n— 1}) is invariant under all unitary operators Ui and thus has to be equal to the whole Hilbert space 
je, i.e. 

n-l 

= 0Ker(XA-A5"a/). (51) 
Therefore, X\ is diagonalizable and can be written in the form 

x, = Y,x^a{u.yp{u^y, (52) 

where P is the projection on the eigensubspace corresponding to the eigenvalue a of Xx and is determined by relations 
of the form 

[/;P(C/;)t = UlP{Ul)^, UJiPiU]:)^ = P ^ 0. (53) 

The equation (|53p applies to an arbitrary pair of unitary operators Ug and Uh in the decomposition of the random 
unitary operation ([T]) and their arbitrary i-th power, i € n. 

The question remains what happens if the set of unitary operators Ui is not irreducible. It is shown in the following 
section that in special cases it may still be possible to decompose the Hilbert space into so-called minimal invariant 
subspaces for which the condition of irreducibility of unitary operators Ui still holds. 



10 



VII. ASYMPTOTIC DYNAMICS OF A TWO-QUBIT CNOT-SYSTEM 



In this section we discuss the asymptotic dynamics of the RUO 

$(p) =PlCipCi + {l~Pl)C2pC2, 



(54) 



which involves two controUed-not (CNOT) operations acting on two qubits. In the computational basis 
{|0, 0), |0, 1), |1, 0), |1, 1)} of the two-qubit Hilbcrt space ^ the action of these CNOTs is defined by 



2\hJ) 



(55) 



with © denoting addition modulo 2. This special RUO of Eq. ([54|l is a hermitian operator and therefore its only 
possible eigenvalue lying within a\i\ are 1 and —1. Let us first find a decomposition of the two-qubit Hilbert space 
into subspaces 1^, i.e. 



(56) 



within each of which the set of unitary operators Ci , C2 acts irreducibly. Constructing such a decomposition is 
equivalent to constructing a decomposition of the finite multiplicative unitary group C generated by Ci and C2. 
The unitary group C is naturally a unitary representation Iq of itself. Therefore, the following considerations are 
immediate consequences of the standard theory of representations of finite groups. 

The unitary group C contains six elements divided into three conjugated classes: Ki = {identity element /}, 
K2 = {C1C2, C2C1}, K3 = {Ci, C2, C1C2C1}. The characters of the representation Iq corresponding to these classes 
are xi = 4, X2 = 1 and X3 — 2. Thus, there are only three inequivalent irreducible representations of the group C, 
say I?^ (fi — {1,2,3}), with dimensions satisfying the relation 



= 6. 



(57) 



Hence, there are two one-dimensional and one two-dimensional inequivalent irreducible representations of the group 
C. The reducible representation Iq can be expressed in terms of irreducible representations as 



(58) 



where are positive or zero integers and fulfil the relation 



3 ^3 
M— 1 -^1—1 



= 5 



(59) 



with g and gi denoting the number group elements and the number of elements of the conjugated class Ki, respectively. 
The only possibility to satisfy the dimensionality of the representation Iq and equation (|59p is the solution: ai = 1 for 
the two-dimensional irreducible representation, a2 — 2 for the one-dimensional irreducible representation, the second 
one-dimensional irreducible representation cannot be involved in the decomposition, i.e. 03 = 0. Two one-dimensional 
representations contained in the irreducible decompositions (|58)) mean that there are just two common eigenvectors for 
the unitary group C and thus common eigenvectors of operators Ci and C2. ^From the definition (j55p it is clear that 
these eigenvectors are ei — |00) and 64 — (|01) + |10) -I- |ll))/\/3. Subsequently, we know that the minimal invariant 
subspaces of operators C'l and C2 are: Vi — span(ei), V2 — span ^62 = :72(l'^"'^^ ^ |10))ie3 — -^(|01) + |10) — 2|11))^ 

and V3 = span(e4). If we denote the restriction of the operator Ci to the subspace Vx as C^^\ in the orthonormal 
basis system {ei}j 



I the operators Ci and C2 correspond to the matrices 
\ 








1 v2 
22" 











\/3 1 

V 0^ 1 / 



C2 - 












a 



(2) 



c. 



(3) 











'\ 







1 
2 


V3 
2 





2 


1 
2 





Vo 








1/ 



(60) 
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Writing the general commutation relations pV.ip in the block structure form we obtain for i € {1, 2} 



(61) 



with the 1x1- matrices X^'^^\ X^'^^\ X^^^\ X'^^^\ with the 2 x 2 matrix X^^^)^ ^j^-jj ^j^g i ^ 2-matrices X^^^) and 
X^^'^\ and with the 2x1 matrices X^^^^ and X'''^'^\ Using equations (pTjl one can check 



^(Ker(X^""-')) C Ker(X^"™0, (Ran(X^""-')) C Ran(X 



(m) , 



(62) 



and thus X'^™"-' is either the zero operator or an invertible operator. Hence, X'^^'^\ X''^'^\ X'''^^\ X'^'^'^'^ are inevitably 
zero matrices. 

Now, assume the case A = 1. A simple evaluation of equation (pT|) leads to the relations X'^'^^^ = a, X'^^'^^ = b, 
^(13) _ ^j^^ X'-^-^'' = d {a,b,c,d G C). The remaining matrix block X'^^^ has to commute with the irreducible 
set of 2 X 2 matrices cl^"* {i € {1,2}) and has to be equal to a multiple of the identity matrix X*^^^^ ^ el [e ^ C). 
The eigenspace of the random unitary operation (|54p corresponding to eigenvalue 1 is five-dimensional and the most 
general eigenvector reads 



/ a c 

e 

e 

Vd 6 



(63) 



The solution of Eq.dni]) with A = -1 yields = X^^^) ^ x^^^^ ^ X^^s) = o. The last matrix block ^(22) jg 

determined by anticommutation relations with the irreducible set of operators Cl {i G {1,2}), i.e. 



C.(2);f(22)^_^(22)^(2)^(22)^ 

^From the discussion in section IVll and the Eq. (|5^ follows that X^^^^ is either the zero operator or 



X(22) ^/(p_Cf^PCf)) 

with the projection P being determined by the equation 

(^(2)^(2)p ^ p^(2)^^2) 

(2) (2) 

Hence, the projection operator P is diagonal in the eigenbasis of the operator C\ . 
general form of the matrix block x(22) corresponding to eigenvalue —1 reads 



(64) 



(65) 



^(22) ^ 



/ 

-/ 



(66) 

Using Eq.(p5|) the most 
(67) 



Thus, the eigenspace of the random unitary operation (j54p corresponding to eigenvalue —1 is one-dimensional and 
the general eigenvector reads 



/O 
0/0 

0-/00 

Vo 



(68) 



Therefore, in the computational basis the attractor space is spanned by the matrices 



/ 1 





Vo 



Xi. 



1 

71 



/O 

11 

10 1 

\0 1 1 



Xi.. 



XlA — 





10 

10 

10 



/O 
10 
10 

Vo 1 



^0 1 


1 
































^,0 





0/ 


























-1 


1 






1 











-1 


1 






(69) 
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In this notation the first index refers to the eigenvalues of the RUO ([54]) and the second index runs through the basis 
states of the corresponding eigenspaces. 

Finally, consider the most general two-qubit input density matrix 



P(0) 



/ an 

''is 

ll4 



V 



ai2 ai3 ai4 

022 023 0,24 

0-23 O33 034 

24 "^34 "^44 



a 



Thus, theorem IIV.3I implies that the asymptotic dynamics under the RUO (E 
determined by the relations 



is periodic with period two and is 



lim p(2n) 

n — >+oo 





a 


c 


c 


c 


c* 


b 


d 


d* 




c* 


d* 


b 


d 


v 


c* 


d 


d* 


b 



lim p{2n+l) = 

n— »+oo 



a 


c 


c 


c 


c* 


b 


d* 


d 


c* 


d 


b 


d* 


c* 


d* 


d 


b 



(70) 



with a = flu, b = ^{022 + 0,33 + 044), c = ^(ai2 + 013 + 014), d = ^(023 + 034 + 024). This two-qubit CNOT network 
is one of the simplest examples of a network allowing for oscillatory asymptotic dynamics. 



VIII. CONCLUSIONS AND OUTLOOK 



We studied general properties of random unitary operations and presented several theorems allowing to determine 
the asymptotic long time dynamics. Thereby, a central result is the structure theorem which states that the asymptotic 
states are located completely inside the vector space spanned by a typically small set of attractors. The form of these 
asymptotic quantum states depends on this attractor space and on the choice of the initial state but is independent of 
the actual values of the probabilities with which the unitary transformations are applied. However, these probabilities 
affect the rate of the convergence towards the asymptotic quantum state. 

It should be stressed that the asymptotic dynamics need not result in a stationary state. Thus, in contrast to 
thermalization the asymptotic dynamics might also be periodic as illustrated by the example of two qubits interacting 
by random C-NOT operations. Even an aperiodic non-stationary asymptotic dynamics is possible. 

The obtained results rise several additional questions. First of all, it is not yet clear what determines the convergence 
rate of a quantum system towards its asymptotic dynamics. Numerical studies suggest that in many cases this 
convergence has an exponential character which depends on the probabilities with which the unitary operations are 
applied. Preliminary results also suggest that at least in the case of many-qubit networks involving controUed-not 
operations the topology of the network is related to the set of attractors. 

Finally, it should also be mentioned that our results might have applications for quantum operations which involve 
an averaging procedure over a group, such as twirling operations. Our results might allow to choose efficiently the 
minimal set of unitary transforms leading to a particular asymptotic state. In addition, we expect that the theory 
presented might also contribute to other related problems concerning the determination of eigenvectors of random 
unitary mapsi2. or their application in purification protocols/ii 
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APPENDIX A: JORDAN CANONICAL FORM 



Let us recall the definition and properties of the Jordan canonical form of square matrices. Consider a complex 
square matrix A = {Aij)^^^-^ of size n x n {A E C"^"). It is similar to a block diagonal matrix 



\ 



A = 



J2 



(Al) 
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in which each Jordan block Ji (i € p := {1, is given by 



J. = 



/A, 1 

A. 

V 



1 

A. J 



(A2) 



Thus, there is an invertible matrix P G C">^" such that A — PJP~^ or equivalently there is a Jordan basis Xi^a G C" 
{i G p, a G {1,2,..., dim( Ji)}) in which the hnear map corresponding to the matrix A has the diagonal form (IA1[) . In 
general, this basis is non-orthogonal and the vectors Xi^a {a G {1,2,..., dim( J^)}) form the basis of the Jordan block 
Ji which corresponds to the eigenvalue Xi of the matrix A. 

The geometric multiplicity of the eigenvalue Xi is the number of Jordan blocks corresponding to Xi and the sum 
of the sizes of all Jordan blocks corresponding to an eigenvalue Xi is its algebraic multiplicity. Therefore, the matrix 
A is diagonalizable if and only if all Jordan blocks are one dimensional. In all other cases any Jordan block, say Ji, 
with dimension s > 1 gives rise to a Jordan chain. This means that there is a so-called lead vector or generator, say 
Xi tijm(ji), which is a generalized eigenvector, i.e. {A — XilYxi^^^^i^j.) =0. The vector i = — Ai/)*^^a;i ciim(Ji) is an 
eigenvector corresponding to the eigenvalue A^. In general, the vector Xij is the image of the vector Xij+i under the 
linear map A — Xil. In this sense all vectors Xi a {ct G {2,3,..., dim( J^)}) are generalized eigenvectors of the matrix 
A. 

Therefore, for every square matrix A there exists a basis consisting only of eigenvectors and generalized eigenvectors 
of the matrix A in which the matrix A can be put in Jordan normal form (|A1[) . 



APPENDIX B: DYNAMICAL MATRICES 



Let us summarize the concept of dynamical matrices which is useful to understand problems related to complete 
positivity of maps. We just recall its definition and present a short summary of characteristic properties needed in 
the main body of our text. Detailed proofs are given in Ref. [7], for example. 

Assume that A is an operator acting on a d-dimensional Hilbert space Jifd- Hence Aij G d :— {l,...,d}) are 
its matrix elements with respect to a given orthonormal basis. It is convenient to interpret a, d x d-matrix {A)ij as a 
vector A = (Am^) G Jf^i^ of the length 

A = (An, ^12, . . . , Aid, ^21, ^22, ■ • • , A2d, . . . , Adi,Ad2, ■ • • , Am). (B1) 

One can check that two d x d matrices A and B fulfil 

{A,B) =Tt{A^B} = A*B ^ {A,B). (B2) 

The vector A of the length d^ may be linearly transformed into the vector A — CA hy a, matrix C of size d^ x d^ 
whose matrix elements may be denoted by C^.^,' with k,k = 1, . . . , d^ . In addition, it is also convenient to use a four 
index notation C™J^ with respect to a two index notation of vectors (|B1[) with m,n, ijl,v — 1, . . . , d^. The matrix 
C may represent an operator acting in a composite Hilbert space = J^d "X" J^d- The tensor product of any two 
orthonormal basis systems in both factors provides a basis in so that we obtain 

CTJ' ^ {em^U\C\er,<E>U) (B3) 

with Latin indices referring to the first subsystem, — .y^d, and Greek indices to the second subsystem, — 
The operation of partial trace over the second or first subsystem produces the dx d matrices = Tr^C or 
= TiA C, respectively, i.e. 

d^ 

and C;^, = (B4) 

m— 1 

If C — B, then C^^l^ = AmnBi^^^^. The standard product of three matrices can be rewritten in the foUowing useful 
form 



u=l 



ABC = $B with $ = A C"^. 



(B5) 
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With the help of identity (IB5[) we can rewrite the RUO ([T|) in the form 

m 

Here, the RUO $ is not understood as a map acting on the d x d-dimensional matrix space but as a map acting on 
the vector space of the dimension cP. 

Let $ be a completely positive trace-preserving map mapping an arbitrary d x d density matrix p € B{J^d) of a 
d-dimensional Hilbert space J^d on a density matrix p e B{Jffd), i-e. 

P = or pI^ = J2 '^Z'Pn.- (B7) 

n,!/— 1 

The meaning of complete positivity becomes rather transparent if we reshuffle $ and define the dynamical matrix D$ 

{D^yp = <i>r • (B8) 
The dynamical matrix £)$ uniquely determines the map $ and has the following properties 

(i) p ^{p'y ^ d^ = di 

(ii) Trp' = 1 ^ Tr^Z?* ^ I mgN 

(iii) ^{I) = I (unital) ^ Trs i:** = / 

(iv) $ is CP map £)$ is positive. 
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